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Abstract 

A  phyllosilicate  is  a  sheet  of  silicate  tetrahedra  bound  by  basal  oxygens.  A  phyl- 
losilicate  automaton  is  a  regular  network  of  finite  state  machines  —  silicon  nodes 
and  oxygen  nodes  —  which  mimics  structure  of  the  phyllosilicate.  A  node  takes 
states  0  and  1.  Each  node  updates  its  state  in  discrete  time  depending  on  a  sum 
of  states  of  its  three  (silicon)  or  six  (oxygen)  neighbours.  Phyllosilicate  automata 
exhibit  localizations  attributed  to  Conway’s  Game  of  Life:  gliders,  oscillators,  still 
lifes,  and  a  glider  gun.  Configurations  and  behaviour  of  typical  localizations,  and 
interactions  between  the  localizations  are  illustrated. 
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1  Introduction 


Phyllosilicates  are  parallel  sheets  of  silicate  tetrahedra,  they  are  widely  present 
in  nature  and  typically  found  in  clay-related  minerals  on  the  Earth  surface 
(I2ll6lf7l23] .  Phyllosilicates  are  used  in  development  of  nano-materials,  nano¬ 
wires  and  patterned  surfaces  for  nano-biological  interfaces  phheh.  They  are 
also  employed,  in  a  form  of  cation-exchanged  sheet  silicates,  as  catalysts  in 
chemical  reactions  |2],  e.g.  nickel  phyllosilicate  catalysts 


Silicate  lattice  exhibit  a  range  of  defects,  or  localised  imperfections.  They  are 
lattice  vacancies,  or  vacancy- impurity  pairs,  interstitial  defect,  and  substitu¬ 
tional  defects  [25H221.  Also  recombination  induced  defects  might  form  when 
a  strong  electrical  pulse  is  applied  to  the  lattice  [28j.  The  defects  can  mi¬ 
grate  through  the  lattice,  especially  in  presence  of  ionising  radiation  I3IB5D], 
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Fig.  1.  Phillosilicate  automata.  Automata  representing  silicon  are  small  discs,  oxy¬ 
gen  automata  are  large  discs.  Unshared  oxygen  atoms  of  the  tetrahedra  are  not 
shown. 

and  interact  with  other  defects,  sometimes  forming  second-generation  defects 
123;  see  also  [22].  In  HE]  we  introduced  phyllosilicate  automata  and  studied 
their  excitable,  threshold  and  intervals  of  excitation  [4],  and  binary-state  [3]. 
By  sampling  of  thousands  of  rules  of  binary  phyllosilicate  automata  [3]  we 
discovered  rules  that  show  (1)  sub-linear  growth  of  disturbances,  and  (2)  ex¬ 
hibit  stationary  and  (3)  travelling  localizations.  As  far  state-transition  rules  of 
phyllosilicate  binary  automata  satisfy  conditions  (1-3),  they  an  be  classed  non- 
rectangular  lattice  analogs  of  Conway’s  Game  of  Life  cellular  automata  [TO], 
We  continue  lines  of  enquiry  into  space-time  dynamics  of  cellular  automata 
on  non-orthogonal  and  aperiodic  lattices,  including  triangular  tessellations 
and  Penrose  tilings  mm.  We  introduce  and  study  an  automaton  network 
-  phyllosilicate  automata  —  where  connections  between  finite  automata  are 
inspired  by  simplified  structure  of  silicate  sheets,  lattices  of  connected  tetra¬ 
hedra,  and  investigate  the  space-time  dynamics  of  compact  clusters  of  non¬ 
resting  nodes  on  these  automaton  networks. 

The  paper  is  structured  as  follows.  In  Sect.[2]we  define  phyllosilicate  automata 
and  Life-like  rules.  Three  types  of  gliders  are  detailed  in  Sect.  [3]  Section  [| 
introduces  stationary  oscillating  localizations.  Stationary  still  localizations  are 
studied  in  Sect.  [5] 


2  Phillosilicate  automata 


Phyllosilicates  are  sheets  of  coordinated  (SiCL)4  tetrahedra  units.  Each  tetra¬ 
hedron  shares  its  three  corner  basal  oxygens  of  neighbouring  tetrahedra  |20ll6l7| 
(Fig.  0-  Vacant  oxygens  of  all  tetrahedra  point  outside  the  lattice,  away  of 
tetrahedra;  these  apical  oxygens  are  not  taken  into  account  in  our  automaton 
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models. 


A  phyllosilicate  automaton  A  is  a  two-dimensional  regular  network  (Fig.  [Tj) 
finite  state  machines,  automata;  we  call  them  nodes.  There  are  two  types  of 
nodes  in  the  network:  silicon  nodes  s  (centre  vertex  of  each  tetrahedron)  and 
oxygen  nodes  o  (corner  vertices  of  each  tetrahedron).  A  silicon  node  has  three 
neighbours.  An  oxygen  node  has  six  neighbours  (Fig.  [Tj) . 

The  nodes  take  two  states  0  (resting  state)  and  1  (non-resting  state)  and 
update  their  states  s 1  and  ol  simultaneously  and  in  discrete  time  t  depending 
on  states  of  their  neighbours.  Let  A  (a)  be  a  number  of  1-state  neighbours  of 
node  a,  0  <  cd(s)  <  3  and  0  <  cd(o)  <  6.  Let  o  and  s  be  instances  of  oxygen 
and  silicon  nodes,  and  o  and  s  be  binary  strings  of  seven  and  four  symbols 
each;  a[i]  is  a  symbol  at  ith  position  of  string  a.  A  node  a,  a  =  o  or  s,  updates 
its  state  by  the  rule: 

at+1  =  aat[a(aY},  (1) 

where  a  =  s  or  o  and  a  =  s  or  o.  When  referring  to  any  particular  rule 
(s0,Si,o0,Oi),  where  string  a,  determines  next  state  of  automaton  being  in 
state  i ,  we  will  be  using  decimal  representations  of  strings  as 

lZ(dec( s0),  dec(si),  dec( o0),  dec( si)) 

.  We  assume  state  0  is  a  quiescent,  or  resting,  state:  a  node  in  state  0  will 
not  take  state  1  if  all  its  neighbours  are  0.  We  will  be  calling  nodes  in  state  0 
resting  nodes,  and  nodes  in  state  1  non-resting  nodes. 

In  our  previous  paper  [3]  we  studied  over  100K  rules,  randomly  sampled,  of 
binary-state  phyllosilicate  automata.  We  uncovered  principle  morphologies  of 
patterns  emerging  and  classified  rules  based  on  shape  of  patterns  generated  by 
the  rules  and  internal  morphologies  of  the  patterns,  and  approximated  distri¬ 
bution  of  functions  on  density  of  patterns  and  activity  of  cell-state  transitions 
and  maps  of  oxygen  and  silicon  densities.  Using  shapes  of  patterns  as  a  key 
characteristic  we  selected  five  classes:  octagonally  shaped  patterns,  almost  cir¬ 
cularly  shaped,  stationary  patterns  non-propagating  beyond  convex  or  concave 
hulls,  patterns  exhibiting  dendritic  growth,  and  patterns  showing  still  local¬ 
izations  and  oscillations.  Internal  morphology  driven  taxonomy  is  comprised 
of  solid  patterns,  labyrinthine  patterns,  wave-like  patterns,  disordered  patters; 
and  still,  travelling  and  propagating  localizations. 

The  100K  were  mapped  onto  activity  versus  density  space  [3],  where  activity 
reflects  a  number  of  nodes  changing  their  states  each  time  step  and  density  is  a 
number  of  non-quiescent  nodes  in  a  fixed-size  domain.  Exploring  the  ’activity- 
density’  space  in  the  loci  with  high  activity  and  low  density  we  discovered  three 
minimal  —  minimality  in  number  of  ’l’s  in  strings  s  and  o  —  rules  that  support 
gliders,  7765  =  77(5, 2, 16,  65),  7768  =  77(5, 2, 16,  68),  and  7772  =  77(5,  2, 16,  72). 
The  rules  have  the  same  rules  for  state  transitions  of  silicon  nodes,  s0  =  (0101) 
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Table  1 

Parameters  of  localizations.  For  each  localisation  L  we  indicate  its  period  p,  mini¬ 
mum  UJmin  and  maximum  u>max  numbers  of  non-resting  nodes,  minimum  u>Tn  and 
maximum  'u;°lax  numbers  of  non-resting  oxygen  nodes,  and  minimum  u^iri  and  max¬ 
imum  w4iax  numbers  of  non-resting  oxygen  nodes.  Minimal  values  in  min  columns 
and  maximum  values  in  max  columns  are  underlined. 
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and  Si  =  (0010),  and  oxygen  nodes  in  resting  states,  o0  =  (0010000).  Only 
behaviour  of  oxygen  nodes  in  non-resting  state  differs:  Oi  =  (1000001)  in  rule 
7^65,  (1000100)  in  TZqs,  and  (1001000)  in  7^72.  Each  of  thee  three  rules  7^65,  7^68 
and  7^72.  generate  its  own  unique  glider  and  several  oscillators.  Parameters  of 
gliders  and  oscillators  are  shown  in  Tab.  [l]  We  have  also  discovered  one  still 
life,  observed  in  rule  1Z 68  and  a  glider  gun,  found  in  rule  TZ72- 


Gliders,  oscillators  and  still  lifes  are  characteristic  patterns  of  Conway’s  Game 
of  Life  cellular  automata.  Phyllosilicate  analogs  of  the  Game  of  Life,  or  Life¬ 
like  automata,  obey  the  following  rules: 


•  Silicon  node  in  state  0  takes  state  1  only  if  it  has  one  or  three  neighbours 
in  state  1. 

•  Silicon  node  in  state  1  remains  in  state  1  if  it  has  two  neighbours  in  state  1. 

•  Oxygen  node  in  state  0  takes  state  1  only  if  it  has  two  neighbours  in  state 

1. 

•  An  oxygen  node  in  state  1  remains  in  state  1  if  it  has  no  neighbours  in  state 
1  or  six  in  rule  7^65,  four  in  rule  7 Ze8  and  three  in  rule  TZ72  neighbours  are 
in  state  1. 
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Fig.  2.  Glider  Ggf  in  rule  Ties-  Glider  moves  west.  Resting  nodes  are  shown  by  dots. 
Non-resting  silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes 
are  shown  by  large  circles. 

3  Gliders 


Glider  G^  is  observed  in  configurations  produced  by  rule  7^65,  see  examples  of 
gliders  travelling  west  and  south-east  in  Figs.  [2]  and  [3]  Glider  is  theaviest 
amongst  gliders  discovered  (Tab.  [I]).  It  has  longest  period  —  16  time  steps  - 
and  maximum  weight  of  13  non-resting  nodes  (Figs.  II  and  [dj)).  The  glider 
undergoes  thickening  and  lengthening  phases.  In  thickening  phase  the  glider 
is  larger,  in  terms  of  maximum  between  non-resting  nodes,  along  a  normal  to 
its  velocity  vector,  see  e.g.  Figs.  [2^,bfjko.  In  lengthening  phase  glider  is  larger 
along  parallel  to  its  velocity  vector,  see  e.g.  Figs.  H-lh.il.  Ratio  of  silicon  to 
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Fig.  3.  Glider  Gg®  in  rule  77-65.  Glider  moves  south-east.  Resting  nodes  are  shown  by 
dots.  Non-resting  silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen 
nodes  are  shown  by  large  circles. 
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Fig.  4.  Head  on  collision  of  gliders  in  rule  77-65,  zer°  vertical  shift.  Glider  travelling 
east  collides  with  glider  travelling  west.  Resting  nodes  are  shown  by  dots.  Non-rest¬ 
ing  silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are 
shown  by  large  circles. 
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Fig.  5.  Glider  Gg|  in  rule  77-68 •  Glider  moves  down  from  north  straight  to  south. 
Resting  nodes  are  shown  by  dots.  Non-resting  silicon  nodes  are  shown  by  small 
black  discs.  Non-resting  oxygen  nodes  are  shown  by  large  circles. 

oxygen  non- resting  nodes  changes  during  the  glider’s  cycle  as  follows:  |,  2,  |, 

442644444248010  J  2 
2’  5’  5’  4  ’  4’  5’  4i  4  ’  4  ’  6’  5’  2  ’  9’  3  4’ 

Most  collisions  between  gliders  in  rule  7Z65  lead  to  explosions:  an  unlimitedly 
growing  pattern  is  formed  in  the  result  of  gliders  interaction.  Head  on  collision 
with  nil  vertical  or  horizontal  shift  between  gliders  is  a  rare  case  of  collision 
where  gliders  do  not  explode  (Fig.  [4j  but  fuse  into  an  oscillator.  One  glider 
travels  east,  left  non- resting  pattern  in  Fig.  [4]a,  another  glider  travels  west, 
right  non-resting  pattern  in  Fig.Jdjn  The  gliders  come  into  ‘contact’  with  each 
other  (Fig.  [4^3)  and  form  a  transient  pattern  (Fig.  [4jo)  of  4  non-resting  silicon 
nodes  and  10  non-resting  oxygen  nodes.  This  pattern  is  then  transformed  into 
an  oscillator  of  period  six  (Fig.  [4jT— j ) . 

Glider  in  rule  77-68  is  shown  in  Fig.  [5j  The  glider  has  period  12,  and  maximum 
weight  13  nodes  (Tab.  [I]).  Glider  reflects  along  its  south-north  axis  during 
motion,  in  middle  of  its  cycle,  namely  glider  state  at  time  step  t,  t  —  1,  2,  •  ■  • , 
is  vertically  symmetric  to  the  glider  state  at  time  step  t  +  6.  The  glider  does 
not  show  a  pronounced  thickening  or  lengthening  but  rather  breathing  around 
its  centre  of  mass.  Ratios  of  non-resting  silicon  nodes  to  non-resting  oxygen 
nodes  in  the  glider’s  cycle  are  changed  as  follows:  |,  |,  |,  |,  |,  |,  |,  |,  |,  |, 

6  2  7 

3’  5’  6' 

Head  on  collision  between  two  gliders  in  lZe8  without  offset  leads  to  explosion 
of  both  gliders.  For  certain  offsets  one  glider  can  survive  collision.  Such  sce- 
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Fig.  6.  Interaction  of  gliders  Ggg  in  rule  'Rgs-  Glider  moving  north  brushes  with  glider 
moving  south.  Resting  nodes  are  shown  by  dots.  Non-resting  silicon  nodes  are  shown 
by  small  black  discs.  Non-resting  oxygen  nodes  are  shown  by  large  circles. 

nario  is  illustrated  in  Fig.  [6}  A  glider  moving  north  (left  non- resting  pattern 
in  Fig.  |6^)  clips  a  glider  moving  south  (right  non-resting  pattern  in  Fig.  |6^). 
The  glider  moving  north  does  not  suffer  any  damage  and  continues  its  jour¬ 
ney  undisturbed.  The  glider  moving  south  becomes  tangled  (Fig.  [6Jd  d)  and 
explodes  (Fig.  [6^-i  and  j). 

Glider  in  rule  7Z72  is  shown  in  Fig.  [7J  It  is  smallest  amongst  gliders  discovered 
and  has  smallest  period.  It  has  just  four  different  states.  Ratios  of  non-resting 
silicons  to  non- resting  oxygens  in  the  glider’s  cycle  are  f ,  §  and  |. 

We  did  not  find  glider  guns  in  rules  7Z65  and  7Z68,  we  believe  they  exist  though 
and  just  wait  to  be  found  by  someone  patient.  Glider  gun  in  rule  7Z72  is  shown 
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Fig.  7.  Glider  G|2  in  rule  77-72.  The  glider  travels  east.  Resting  nodes  are  shown  by 
dots.  Non-resting  silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen 
nodes  are  shown  by  large  circles. 
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Fig.  8.  Glider  gun  in  rule  77-72.  Resting  nodes  are  shown  by  dots.  Non-resting  silicon 
nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown  by  large 
circles. 

in  Fig.  [8] 


4  Oscillators 


Oscillator  is  a  localised  compact  configuration  of  non-quiescent  states  which 
undergoes  finite  growth  and  modification  but  returns  to  its  original  state  in  a 
finite  number  of  steps.  Oscillators  are  generated  by  all  three  rules,  discussed 
in  the  paper,  and  can  be  easily  found  by  randomly  perturbing  the  automaton 
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Fig.  10.  Oscillator  Og5  in  rule  7^65-  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 
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Fig.  11.  Oscillator  Ogg  in  rule  1Zq§.  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 


lattice.  Here  we  present  a  selection  of  most  common  oscillators.  Some  param¬ 
eters  of  the  oscillators  are  shown  in  Tab.  [lj  Periods  of  oscillators  vary  from 
2,  oscillator  Oj2,  to  12,  oscillators  Og|  and  0?$.  Minimum  weights  vary  from 
2,  oscillator  0)j\  to  6,  oscillator  0|5.  Maximum  weights  vary  from  4,  oscillator 
Of2,  to  26,  oscillator  Ogg. 


Three  typical  oscillators  in  rule  7?.65  are  shown  in  Fig.  |9j[T0]and  11  Oscillator 


Og5  changes  it  configuration  from  two  parallel  chains  of  non-resting  oxygen 
nodes  oriented  north-west  to  south-east  (Fig.  [9^)  and  north-east  to  south¬ 
west  (Fig.  [9]:)  to  sparse  chains  of  non-resting  silicons  with  few  non-resting 
oxygens  incorporated  (Fig.  [^d).  Oscillator  0|5  is  also  oscillator  in  rule  1Z 72. 
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Fig.  12.  Oscillator  Og8  in  rule  TZes-  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 


o  o 

o  o 

o  o 

o  o 

m 

:o  9f  o:  : 

:  .  O  A  O  ; 

:  ^o*0- 

•O.  °  V3* 

of  Oo 
° 

(a)  t 

(b)  t  +  1 

(c)  t  +  2 

(d)  t  +  3 

(e)  t  +  4 

•  oooo  • 
.  oooo  . 

*00* 


(f)  t  +  5 


oooo 
o  oQo  o 


o  o°o  o 
oooo 


9  :  9 

•  o  * 


oo 


oo 


oo 


oo 


(g)  t  +  6  (h)  t  +  7  (i)  t  +  8 


oooo 

o°o 

°o° 

...  ... 

oo 

o&oo 

o  o 

o  o 

•  •  •  •  •  • 

oo 

(k)  t  +  10  (£)  t  +  11  (m)  t  +  12 


oo 

.  :  oo  : 

*  ^  * 

.  . 

■  ■  ■  oo  : 

:  '  OCl  ■ : 

(j)  t  +  9 


Fig.  13.  Oscillator  Og|  in  rule  7^g§.  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 


Oscillator  Og5  switches  its  configuration  between  a  single  chain  of  non-resting 
silicons  (Figs.  [Tojig).  a  single  chain  of  non-resting  oxygens  (Figs.  [Toja),  and 
regular  non-linear  arrangements  of  non-resting  nodes  (Figs.  lCb-f).  Oscillator 
is  one  of  smallest  oscillators  (Fig.  11).  It  has  period  six.  The  oscillator 


configuration  at  time  step  t  is  symmetric  —  along  north-east  to  south-west 
to  its  configuration  at  time  step  t  +  3. 


axis 


Exemplar  oscillators  discovered  in  rule  are  shown  Figs.  [12]  and  13 
cillator  Og8  has  the  same  configuration  as  oscillator  see  Fig.  12 
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Fig.  14.  Oscillator  Oy2  iu  rule  ^72-  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 
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Fig.  15.  Oscillator  Oy2  in  rule  7^72-  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 
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Fig.  16.  Oscillator  0\ \  in  rule  7^.72-  Resting  nodes  are  shown  by  dots.  Non-resting 
silicon  nodes  are  shown  by  small  black  discs.  Non-resting  oxygen  nodes  are  shown 
by  large  circles. 

Fig.  [njicij.  Oscillator  Og|  is  the  largest  oscillator,  amongst  presented  in  the 
paper  (Tab.  [l]).  The  oscillator  has  period  16,  and  consists  of  28  non- resting 
nodes  in  its  heaviest  configuration  (Fig.  [T3]f) . 

Oscillator  0^2  is  the  smallest  one  (Tab.  [TJ.  It  has  just  two  configurations: 
one  consists  of  four  non-resting  silicon  nodes  and  another  of  three  non-resting 
oxygen  nodes  (Fig.  |l4|). 


Oscillator  0} \  (Fig.  16)  is  the  smallest  oscillator  with  longest  period.  Its  period 
is  12  time  steps  and  maximum  weight  is  5  non-resting  nodes,  see  e.g.  Fig.  16  i. 
When  glider  Gf2  collides  to  O] \  there  are  three  possible  outcomes,  depending 
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(b)  t  +  1  (c)  t  +  2  (d)  t  +  3  (e)  t  +  4 


Fig.  17.  Still  localisation,  (a)  Still  localisation  in  rule  77-68;  in  this  rule  the  localisation 
acts  as  eater.  The  localisation  (a)  becomes  an  oscillation  (b-e)  in  rule  77 65.  Resting 
nodes  are  shown  by  dots.  Non-resting  silicon  nodes  are  shown  by  small  black  discs. 
Non-resting  oxygen  nodes  are  shown  by  large  circles. 

on  state  of  the  oscillator  at  the  collision  and  relative  positions  of  the  glider 
and  the  oscillator:  oscillator  and  glider  vanish  in  the  result  of  collision,  both 
localizations  explode,  and  glider  disappear  but  oscillator  remains  intact.  In 
Conway’s  Game  of  Life  an  eater  is  a  still  life  (stationary  localisation  that  does 
not  change  its  configuration)  which  repairs  itself  after  collision  with  a  glider 
or  any  other  travelling  localisation,  while  the  colliding  travelling  localisation 
vanishes.  Oscillator  0^\  is  not  a  still  life  so  it  can  be  classified  as  oscillator- 
eater. 


5  Still  life 


Still  lifes  are  localizations  that  do  not  change  their  state/pattern  or  a  position 
during  automaton  development.  So  far  we  found  still  lifes  only  in  rule  776s:  this 
is  a  ring  of  six  oxygen  and  six  silicon  non-resting  nodes  as  shown  in  Fig.  [T7}.i. 
This  patterns  becomes  an  oscillator  in  rule  77-65  (Fig.  [T?)))  and  is  transformed 
into  a  glider  gun  in  rule  7768  (Fig.  [8]). 


6  Discussion 


We  introduced  a  regular  automaton  network  which  mimics  topology  of  a  sil¬ 
icate  sheet  and  studied  node  state  transition  rules  which  exhibit  behaviour 
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Fig.  18.  Positions  of  Game  of  Life  and  phyllosilicate  automata  rules  on  ‘neigh¬ 
bourhood  occupancy  necessary  for  birth’,  transition  from  state  0  to  state  1,  versus 
‘neighbourhood  occupancy  necessary  for  birth’,  transition  from  state  1  to  state  1. 
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similar  to  Conway’s  Game  of  Life  automata. 


Classical  Conway’s  Game  of  Life  automata  have  rule  B3/S23:  a  dead  cell 
become  alive,  born  if  it  has  three  live  neighbours,  a  living  cell  survives  if  it  has 
two  or  three  live  neighbours.  In  notation  of  Life-like  automata  the  rules  can 
be  written  as  follows.  All  silicon  nodes  update  their  states  by  rule  B23/S23. 
Rules  of  oxygen  nodes  are  B2/S06  (rule  72.6s),  B2/S04  (rule  Ties)  and  B2/S03 
(rule  7^72).  Such  notations  could  be  somewhat  misleading.  In  Conway’s  Game 
of  Life  a  node  (cell)  has  eight  neighbours  while  in  phyllosilicate  automata  a 
silicon  node  has  three  neighbours  and  oxygen  node  six  neighbours.  Thus,  rules 
should  be  normalised  by  dividing  number  of  non-resting  neighbours  necessary 
for  birth  and  survival  by  neighbourhood  size.  Normalised  rules  are  as  follows. 
Normalised  Conway’s  Game  of  Life  rule  is  B|/S||.  Rule  for  silicon  nodes  in 
phyllosilicate  is  B|l/S|.  Rules  for  oxygen  nodes  are  B|/S01  (72.6s),  B|/S0| 
(Ties),  and  B|/S0|  (72.72).  On  the  map  ‘neighbourhood  occupancy  necessary 
for  birth’,  transition  from  state  0  to  state  1,  versus  ‘neighbourhood  occupancy 


necessary  for  birth’  (Fig.  18)  Game  of  Life  rules  reside  in  the  quadrant  with 
below  average  occupancies.  Occupancies  necessary  for  silicon  nodes  to  take 
non-resting  state  and  to  remain  in  the  non-resting  reside  in  the  quadrant  of 
above  average  occupancies.  Rules  for  oxygen  nodes  are  characterised  by  below 


average  occupancies  necessary  for  survival  (Fig.  18).  Distribution  of  rules  on 
this  map  hints  that  conditions  for  survival  are  more  critical  than  conditions  for 
birth.  Most  rules  supporting  gliders  obey  a  neigbourhood  occupancy  between 
0.3  and  0.4  while  neighbourhood  occupancy  for  birth  varies  from  0  to  1. 


We  presented  three  types  of  node-state  transition  rules  of  a  regular  automaton 
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network  which  mimics  topology  of  a  silicate  sheet.  All  three  rules  studied 
support  travelling  and  oscillating  stationary  localizations.  Yet  only  one  rule 
exhibits  still  localizations.  Our  search  for  phyllosilicate  automata  exhibiting 
Conway’s  Game  of  Life  behaviour  was  extensive  yet  not  exhaustive.  Thus 
substantial  chances  remain  that  glider  guns  could  be  found  in  rules  1Z 65  and 
Ties,  and  still  lifes  in  rules  7^65  and  TZ^- 

Most  interactions  between  localizations  lead  to  explosions  of  the  localizations, 
when  unlimitedly  growing  patterns  of  non-resting  states  are  formed.  Just  few 
scenario  of  collisions  lead  to  annihilation  of  both  or  at  least  one  of  colliding 
localizations.  This  may  somehow  limit,  however  not  totally  take  away,  an  ap¬ 
plicability  of  phyllosilicate  automata  in  design  of  collision-based  circuits  [I]. 
Further  studies  in  the  interaction  of  localizations  are  required  to  make  a  def¬ 
inite  conclusion  about  usefulness  of  phyllosilicate  automata  in  the  unconven¬ 
tional  computing  field. 

How  gliders,  oscillators  and  still  lifes  could  be  represented  in  real  phyllosilicate 
lattices? 

We  believe  they  will  be  seen  as  domains  of  localised  defects,  e.g.  recombination- 
induced  defects  formation  [28j  occurred  in  tetrahedra,  and  vacancy-impurity 
pairs  or  isolated  interstitial  silicon  defects  [30] .  The  imperfection,  vacancy, 
and  self-interstitial,  substitutional  and  interstitial  impurities  are  important 
because  they  are  necessary  for  setting  up  travelling  localizations  of  e.g.  a  sub¬ 
stitutional  impurity  atom  through  the  lattice.  The  travelling  impurities  cause 
local  deformations  of  the  phyllosilicate  lattice,  e.g.  via  recombination  of  vacan¬ 
cies  and  carriers,  these  deformations  produce  waves.  Thus,  for  example,  state 
T’  of  a  node  can  represent  a  point  interstitial  defect,  where  the  host  atom 
moves  to  a  non- lattice  position  [25].  The  point  defects  can  form  travelling  lo¬ 
calizations,  phenomenologically  similar  to  gliders.  Such  localizations  of  defects 
may  not  be  stable,  even  at  low  temperature,  and  disappear.  However,  dura¬ 
tion  of  their  life  time  could  be  enough  to  execute  one  or  more  collision-based 
logical  gates. 

A  typical  scenario  of  generating  a  glider  on  a  phyllosilicate  lattice  would  be 
the  one  detailed  in  ra.  A  locally  irradiation  of  the  lattice  with  1-3  MeV 
electrons,  gamma  rays,  or  fast  neutrons  |T3]  would  move  silicon  atoms  into 
metastable  interstitial  states.  Thus  a  negatively  charged  travelling  vacancy  is 
formed.  Direction  of  this  localisation  travelling  could  be  controlled  by  inducing 
lattice  vibrations  and  affecting  electrons  transitions.  Obviously,  to  obtain  an 
experimental  confirmation  would  be  a  very  difficult  task. 

There  is  at  least  some  analogy  between  gliders-supporting  automata  rules 
and  conditions  for  emergence  of  a  travelling  vacancy.  In  all  glider-supporting 
phyllosilicate  automata  studied  a  silicon  node  excites  and  stays  excited  if 


15 


it  has  two  or  three  excited  neighbours.  An  oxygen  is  excited  if  it  has  two 
excited  neighbours.  As  hinted  in  [T3j  negatively  charged  travelling  vacancies 
can  migrate  only  if  at  least  two  inter-atomic  bonds  are  broken. 

When  complemented  by  molecular-dynamics  method  [ITj  or  reaction-diffusion 
computational  approach  ra  of  simulating  defects  propagation,  our  approach 
can  be  useful  in  designing  close  to  reality  computing  circuits  based  on  travel¬ 
ling  defects  and  clusters  of  defects.  This  will  be  a  topic  of  future  studies  to  find 
what  exact  structure  of  defect,  localisation  or  dislocation  ISEIEIj  corresponds 
corresponds  to  gliders  and  still  lifes. 
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